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Abstract 

In this paper we propose five different modifications of the hard sphere potential for the 
modeling a short-range repulsion and the calculation of thermodynamic and transport properties 
of liquid 4 He. We calculate the potential energy, the total energy, and the sound velocity at 
T = K. It is shown that three of the proposed potentials give a satisfactory description of 
these properties. 
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1 Introduction 

In this paper we calculate some properties of the helium system using the modifications of helium- 
helium interatomic potential obtained in |]. 

The method of receiving the potential is based on the collective variables formalism as described 
in our previous work [|1|]. 

Numerical methods could be grouped in two large parts: 1) calculations with adjustable param- 
eters using experimental data on thermodynamic and transport properties [Q, |3|, ^, [J, 2) ab initio 
quantum- mechanical computations (mainly Monte Carlo) || [7j. A good comparison between various 
potentials, calculated properties of helium, and the experimental data might be found in ||. 
>• ■ We analyze our potentials due to the form of the potential curve and such properties as the 
potential energy, the total energy, and the sound velocity, all calculated at T — K, i. e., for the 
^ \ ground state. 

Concerning the form of the potential curve, one can notice that the short-range part of the 
potential is not very significant starting from distances about 2 A and less. This statement becomes 
clear if one looks at the pair distribution function F 2 (R) of helium ||: F 2 {R < 2.212 A) = 0. 

On the other hand, the results for the potential from PJ lead to the value — 7.8 K for the total 



energy of the system of helium atoms at K, while the experiment gives — 7.17 K [|Kj, we consider it 
as a good agreement. But we obtained the value of the potential energy —15.5 K, and computer sim- 
ulations |§ predict it in the range from about —20 to —22 K (there is no experimental measurements 
available for this quantity). This difference may not seem to be quite satisfactory if one calculates 
quantities connected with the potential energy only. 
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Thus, we tried to correct our potential for it could reflect the short-range interactions more 
precisely unlike that obtained in fll]. For this purpose we used some modification of the hard sphere 
potential. 

The paper is organized as follows. The calculating procedure is described in the Section 2. 
Different forms for the short-range repulsion are presented in the Section 3, the respective potential 
curves are also given there. In the Section 4 we adduce the computation of the energy and the sound 
velocity of 4 He in the ground state. 

2 Calculating procedure 

The following relation was used for the calculation of the potential Fourier image v q . see []fl |Tj]: 

+ q 2a 2(q)-q 2 a^(q)-^^k 2 a 4 (k,-k,q, -q) (1) 

+ Tf ^ k(_k ~ q ) a ^ q ' k ' ~ k ~ q ) = °' 

N is the total number of particles, V is the volume of the system, in the thermodynamic limit we 
have N/V = const =p = 0.02185 A~ 3 ; m is the mass of 4 He atom; the quantities 03, are expressed 
through a 2 , see JTJ, O. 

We are going to solve this problem in the approximation of "two sums over the wave vector". 
Such calculations were made in Q with the expression for a 2 in the form of 

S(q) being complex functional of a 2 containing summation (integration) over the wave vector Jl|, 11 . 

The equation ([[]) for v q is solved using iterative procedure. We have used the corrected expression 
for the zeroth approximation of the quantity a 2 ||12||: 



while in JTJ we used 

aS(0 = i (l - ; (4) 

Sg stands for the short-range structure factor, i. e. it corresponds to the short-range repulsive 
interactions. The question obtaining this structure factor S® is elucidated in the next section. 

3 Short-range potential 
3.1 Hard sphere repulsion 

In this paper we first consider a short-range potential in the form of hard sphere interaction, 
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a is the hard sphere diameter. 

This leads to S q = S^ s , the structure factor of the hard sphere system. It was calculated from 
the pair distribution function of the hard sphere system in which we neglected all terms with orders 
of density higher than 2. This approximation will be called here as /^-approximation or HS-p 2 . 

It is known that the Percus-Yevick approximation (PY) gives quite good description in the case 
of short-range interactions. It has also another advantage: there is an exact solution of the Ornstein- 
Zernike equation for the hard sphere system in PY ||13|| . We denote the structure factor obtained in 
this fashion as S^ Y . 

The quantities S^ s — 1 and S^ Y — 1 show a weak damping oc 1/q 2 at the large values of the 
wave vector. It leads to the complication of the calculation procedure since one should extend the 
integration limits in the Fourier transformation of v q . More precisely speaking, the damping of v q 
was weak too and we used the extrapolation at large q to obtain correct behaviour of <&(R) at small 
distances. 

We have calculated the Fourier image v q and the potential $(-R) using the hard sphere structure 
factor both in HS-p 2 and PY approximations. As one can see from the Fig. [1], the long-range 
behaviour of this potential is not satisfactory unlike the short-range one. 

One can see from Fig. |l] the extra hill on the potential curve at approximately 6 A. We will try 
to find out the reasons of its annearinEf. 
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Figure 1: The potential $(-R) with short-range interactions in hard sphere approximation (HS-p 2 
solid line, PY — dotted line) compared with results of our previous work Q (dashed line). 



3.2 Softened repulsion 

The true repulsive part of the potential should be a bit "softened". It is easier to formulate this 
"softening" in the terms of Meyer function f(R): 

f(R) = e -0*(«) - 1 (6) 

that stays finite while $(-R) can have infinities, see (|5]). 
For the hard sphere system f(R) reads 

-1, R < a m 
0, R>a ' { ) 

In this terms, the softening means that the straight step of the Meyer function of hard sphere 
must be slightly "rounded" . Here we consider the variants of this "rounding" . 



f w (R) = 6{R-a)-l = 
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Let us suppose the power dependence in $°(i?): oc 1/ R n . Then the kinetic energy operator 

in the Schrodinger equation gives: 

Thus, the power 2n + 2 must equal the power of inverse distance in the repulsive part of real potential. 
Assuming it to be the Lennard-Jones potential we obtain n = 5. 

Another way is to round the corners in the Meyer function via direct accepting its form. For 
instance, like in the previous case: 

f(R) = e~^ n - 1. (9) 

The value of n should be taken large (larger than 5) in order to resemble the hard sphere step. 

The first case (Lennard-Jones repulsion, n = 5) will be called "soft spheres" (SS), and the second 
one (directly accepted repulsion) will be denoted as "almost hard spheres" (AHS). 

Let us choose the following values for AHS quantities in @: 

A = 2.1 A, n = 12. (10) 

The Meyer function then becomes quite close to the hard sphere step, see Fig. That is the reason 
why we call it "almost hard spheres" . 
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Figure 2: Meyer functions of hard spheres (solid line) and almost hard spheres (dashed line). 

We have calculated the structure factor of the softened repulsive part using the PY approximation. 
Since we do not need big accuracy, only the first iteration in the Ornstein-Zernike equation was used. 
Our results show that the potential curve does not differ much from that of our previous work [Ij], 
see Fig. |3|. 

In spite of large difference in 1/R powers (5 and 12), the results do not yield to much different 
behaviour in the region of small R: it still resembles that of 0. From this point of view, the reason 
must be in principal difference between soft and hard spheres, namely in the analyticity or non- 
analyticity of Meyer function that leads to different types of damping at large q for the structure 
factor and, result, for the Fourier image v q . 

Another fact concerns the region of potential well. One can see that the increasing of 1/R 
power from 5 to 12 moves the depth of the well closer to that of obtained from the hard sphere 
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Figure 3: The potential $(-R) with short-range interactions in SS (dashed line) and AHS (solid line) 
approximation compared with results of our previous work [I]] (dotted line) 

approximation. On the other hand, similar value was obtained in |]J and it appears not to differ 
much from the results of computer simulations || [7| and the computations with adjustable parameters 
0, 0, [|, |5|. Those results might be found in Table 1. 



Potential 


$(0), K 


Well 


Extra Hill 


Ref. g 

HS - p 2 
HS - PY 
AHS 
SS 

Glued 


~ 240 

oo 

oo 
~ 240 
~ 230 

oo 


-10.59 K at 3.34 A 
-10.86 K at 3.32 A 
-11.55 K at 3.27 A 
-11.76 K at 3.34 A 
-15.95 K at 3.32 A 
-11.72 K at 3.33 A 


0.26 K at 6.41 A 
1.99 K at 4.84 A 
3.64 K at 4.94 A 
0.23 K at 6.31 A 
0.27 K at 6.44 A 
0.23 K at 6.31 A 


Ref. [f 
Ref. [ 


3 


2 • 10 6 


-11.01 K at 3.0 A 
-11.04 K at 2.99 A 


no hill 
no hill 



Table 1. Parameters of the potential curve for different types of potential. 

The extra-hill was obtained on our potential curves at the interatomic distances of ~ 5 — 7 A. We 
consider it as a consequence of solving the problem of interaction between helium atoms via quantum 
equations. But the large values of hill maxima obtained in the case of hard-sphere interaction should 
be treated as the result of Meyer function non-analyticity. We claim it since the correspondent values 
in the case of soft spheres and almost hard spheres are commensurable with the value from [T| . 

The region of small R values is probably the most complicated one. It seems that one should 
accept the hard-sphere approximation in order to obtain a good description here. We tried to join 
the advantages of the PY short-range part with those of the AHS long-range part by just "gluing" 
this two data sets at the point of 3.25 A. The left side corresponds to PY; the right side - to AHS. 

One can notice that there is a constant positive difference between the minimum position (the 
equilibrium interatomic distance) of the potentials discussed in this work and of those calculated by 
Aziz et al. We consider it as a consequence of the effective nature of our potentials. Let us consider 
the model of hard spheres with diameter a in order to illustrate this statement. The equilibrium 
distance between two atoms in the system of two atoms is obviously a. But if one considers the 
system of three particles than the average distance slightly increases, namely by 5 per cent. Thus, 
the difference in approximately 10 per cent in our case might be at least partially considered as a 
consequence of such a geometric effect. 
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4 Energy and sound velocity 



One can calculate the potential energy $ of the system at K using the relation [14 

' 1 



*="T + S 2 



- 1 



'1 + 



2pv q 



2^2 



hfq_ 
2m 



(11) 
(12) 



The total energy E reads: 



E = $ + 



q^O 



(13) 



The quantity pv§ in these equation must be substituted with mc 2 where c is the sound velocity 
in 4 He. The value mc 2 = pv§ corresponds to the RPA and we want to have the next, postRPA 



approximation. Thus, the corrected expression reads (12 



1 



mc 



a q - 1) 



OLa 



(14) 



We give the results discussed above in the Table 1 for the comparison. The left hand values of 
sound velocity correspond to the RPA (mc 2 = pz/ ), while the right-hand values reflect the correc- 
tion pi. 



Potential 


$/JV, K 


E/N, K 


c, m/s 


Ref. 
HS-p 2 a ) 
AHS 
SS 

Glued 


-15.5 
+18.1 
-21.6 
-29.1 
-16.7 


-7.8 
+38.8 
-12.0 
-18.6 

-8.0 


292 vs 238 
407 vs 524 
297 vs 240 
254 vs 182 
324 vs 262 


Ref. | b ) 


-21.7 


-7.1 


241 


Experiment 




-7.17 c ) 


240.3 d ) 



a ) We do not give data for the potential obtained using S PY since in this case one cannot obtain a "good" a q : value 
of 1 + 2pv q /e q is negative for some values of q. 

b ) Calculation using potential 0]. 

c ) Ref. 0. 

d ) This value corresponds to the sound velocity at saturated vapour pressure for T = 0.1 K Ref. [l6| gives 
237.8 m/s at T — 0.8 K. These values are not strictly experimental but it is possible to consider them as the averaged 
experimental ones since they are calculated using the well-analyzed equation of state. 



Table 2. Energy and sound velocity with different potentials. 



As might be seen from the Table 2, the HS model for the short-range repulsion yields to incorrect 
values for energy and the sound velocity. Unfortunately, we did not find the universal potential sat- 
isfying the requirements to describe the total energy and the potential energy with similar accuracy. 
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The best fit for both these quantities is found in the case of "glued" potential but the method of 
its obtaining looks very artificially. One can see that the best fit of the potential energy is found 
in the AHS approximation for the short-range repulsion. On the other hand, previously obtained 
potential |l| remains the best of those considered here due to the values of the total energy and the 
sound velocity. 

Therefore, three potential will be used in our further works: the potential obtained in the 
potential with AHS short-range repulsion, and the "glued" potential. The comprehensive analysis of 
thermodynamic and transport properties that might be calculated and compared with the experiment 
or the computer simulations will judge the most suitable model potential. 
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